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We investigate the spectroscopy and decays of the charm-strange quarkonium system in a potential 
model consisting of a relativistic kinetic energy term, a linear confining term including its scalar 
and vector relativistic corrections and the complete perturbative one-loop quantum chromodynamic 

■ short distance potential. The masses and wave functions of the various states are obtained using 
| a variational technique, which are then used in a perturbative treatment of the potential to find 

£\j , the mass spectrum of the cs system and radiative decay widths. Our results compare well with the 

, ^ available data for the spectrum of D s states. We include a discussion of the effect of mixing and an 

' investigation of the Lorentz nature of the confining potential. 

s" 

■ Recently we reported on a study of the charmonium and upsilon systems in a semi-relativistic model which 
' includes all v 2 /c 2 and one- loop QCD corrections for the interaction of a quark and antiquark of equal mass[l]. 

This semi-relativistic potential model successfully describes the spectra and leptonic and radiative decays of 
those systems. We have now extended this modelling approach to systems in which the quark and antiquark 
have different masses. 

Interest in the modelling of light-heavy quarkonia is over 25 years old [2] . A variety of modelling approaches 
have been employed with varying success [3-8] . Renewed and continuing interest in the modelling of cs quarkonia 
Jy_j is fueled by, in particular, the recent discovery of the 2 3 Si state [9] as well as ongoing efforts to determine the 

, masses and decays of the D s mesons [10]. 

, We have revised and extended the approach of our earlier papers in order to investigate the spectroscopy and 

' decays of the D s system, as well as to discuss other questions of modelling interest. In addition, we investigate 
, the the scalar/ vector mixture of the phcnomcnological confining potential 
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INTRODUCTION 



In the next Section, we describe the potential model in some detail. This is followed, in Section 3, by an 
outline of our calculational approach. In Section 4, we present our results for the D s system, and then give some 
conclusions in Section 5. The conventions we use for our treatment of the mixing of the J = 1 p— states are given 
in the Appendix. 



SEMI-RELATIVISTIC MODEL 



In our analysis, we use a semi-relativistic Hamiltonian of the general form 



H = y p 2 + m\ + y p 2 + m 2 + Ar — 



l-^ + ^(33-2n / )(ln(/ir)+ 7 ^ 



+ V L + V S (1) 

= H Q + V L + Vs, (2) 



* Electronic address : sr adf ordabr o ckport . edu 
^Electronic address: repkoOpa.msu.edu 
-^Electronic address: mike.saelimagmail.com 



2 



where mi and mi are the quark masses, fi is the renormalization scale, 77/ is the effective number of light quark 
flavors and 7^ is Euler's constant. Vl contains the v 2 /c 2 corrections to the linear confining potential 
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where A is the linear coupling strength. The first line in Eq.(3) is the contribution from scalar exchange while 
the second and third lines give the contribution from vector exchange, with fy representing the fraction of vector 
exchange in the interaction. The short distance potential is [2] 



with 



V s = V HF + V LS + V t + Vsi + Vmix , 



(4) 
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(5d) 



(5e) 



We have chosen H Q such that it contains the relativistic kinetic energy and the leading order spin-independent 
portions of the long-range confining potential and the one-loop QCD short-range potential. It is important 
to recall that the potential given by Eq. (4) does not reduce to the potential in Rcf.l, due to the presence of 
annihilation terms in the equal-mass quark-antiquark potential. It should also be noted that in calculating the 
matrix elements of the 5(r) terms in Eqs. (5a) and (5d), we 'soften' their singularity by adopting the quasistatic 
approximation of Ref. [11], which leads to the replacement 

2 

Tfi 

7rr 

where m is the quark mass. This softening helps the stability of the eigenvalue calculation. 
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3. CALCULATIONAL APPROACH 



The cs mass spectrum and corresponding wave functions arc obtained using the variational approach described 
in Ref.l. In this approach, we expand the wave functions as 

™ / r \ k +? 

Wutf) = E c * (r) e ~ r/R yrM , (7) 

fc=0 

where 3^? s (^) denotes the orbital-spin wave function for a specific total angular momentum j, orbital angular 
momentum £ and total spin s. The CVs are determined by minimizing 

with respect to variations in these coefficients. This procedure results in a linear eigenvalue equation for the 
Cfc's and the energies, and is equivalent to solving the Schrodinger equation. The wave functions corresponding 
to different eigenvalues are orthogonal and the fc th eigenvalue Afe is an upper bound on the exact energy E].. 
For n = 14, the lowest four eigenvalues for any £ are stable to one part in 10 6 . We performed a perturbative 
calculation, using Hq as the unperturbed Hamiltonian and all other terms treated as first-order perturbations. 

An optimal set of the parameters a = (ai, a%, • • • , a n ) were found by minimizing the \ 2 function 

2 _ (Ccxpi - Oth{a)t) 2 , s 

i=l 1 

where the Oi denote the experimental and theoretical values of some quarkonium observable and the Oi are the 
associated errors. In this work, the O e xpi consist of a subset of the measured D s masses. For the masses, the 
<7j are taken to be the actual experimental error and a common intrinsic theoretical error added in quadrature. 
The latter error reflects the theory uncertainty associated with omitting corrections beyond one-loop and is 
estimated by requiring the x 2 /degree of freedom to be approximately unity. Typically, this error is a few McV. 
The minimization of \ 2 with respect to variations of the parameters a is accomplished using the search program 
STEPIT [12]. The choices of as and m c were kept consistent with the results of running these parameters from 
the charmonium scale of Ref. [1] by introducing a Gaussian prior in the \ 2 function for each one. For additional 
discussion of calculational details, see Appendix A of Ref.l. 



4. RESULTS 



We summarize our results in the following tables. The parameters resulting from our fit are given in Table I. 



A (GeV 2 ) 


0.115 


as 


0.391 


m c (GeV) 


1.66 


m s (GeV) 


0.346 


H (GeV) 


1.19 


fv 


0.00 



TABLE I: Fitted Parameters for the cs system 

The results for our determination of the D s levels are shown in Table II [13]. Overall our fit to the spectrum 
is quite good. 

As is usual in potential model treatments [5, 7, 14-16], the radiative widths were calculated in the dipolc 
approximation. We obtained the E\ matrix elements by using the variational radial wave functions to construct 
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m C s (MeV) 


Model 


Expt 


D s 


1968.49 


1968.49 ±0.34 


D* a 


2112.3 


2112.3 ±0.5 


D s0 (2317) 


2317.8 


2317.8 ±0.6 


D s i(2460) 


2459.6 


2459.6 ± 0.6 


Dsi(2536) 


2535.2 


2535.35 ±0.34 


D s2 (2573) 


2572.9 


2572.6 ±0.9 


D 3 (2486) 


2485.8 




D s *(2637) 


2637.2 


2690. ± 7 



TABLE II: Results for the D„ spectrum are shown. The fit uses all the indicated states of the D a system except for the 
n = 2 D 3 and D* . 



initial and final state wave functions with the appropriate angular dependence and explicitly performing the 
angular integration. Our results are 



4a / q 1 m 2 - g 2 "H 
9 \ mi ± m,2 



w 3 K/|r|i)| 



2 E f 

Mi 



1 for 3 Pj 

1 for 1 P 1 

(2J + l)/3 for 3 S 1 

3 for 1 S () 



-> 3 Si 



(10) 



for Ei transitions. Here, u> is the photon energy, q\ and q 2 are the quark charges in units of the proton charge, 
Ef is the energy of the final quarkonium state, Mi is the mass of the initial quarkonium state, and mi and m 2 
are the quark masses. 

We also take into account the mixing between the 1 P\ and 3 Pi eigenstates of the cs Hamiltonian due to the 
L-(S\ — S2) terms in Eqs. (3) and (5e) of the perturbative potential. This mixing yields the two J = 1 states 
D s 1 and D' s 1 . They are, explicitly, 

|D a i(2460)> = sin(6»)| 3 Pi}±cos(6')| 1 Pi), (11a) 



| J D; sl (2536)) = cos(0)| 3 Pi) -sin^Pi), 



where 



tan(0) = - 



V 31 



E. 



(lib) 



(12) 



with V31 denoting the expectation value of the mixing terms and E+ denoting the larger of the two eigenvalues 
of the mixing matrix. Note that, because of the l/m| behavior of the these terms, the mixing is quite sensitive 
to the strange quark mass. The conventions used in parameterizing the mixing arc given in the Appendix. 

For Mi transitions, a parallel calculation, using the fact that the singlet and triplet s— states have the same 
wave functions in the perturbative treatment, yields 



4a , / <?i q 2 \ 2 E f 
^-T^ v2^~2^J M (13) 

for the D* s — > D s ± 7. In the case of the p— state magnetic transitions D s 1 — > D s + 7 and D' sl — > D s + 7, both 
the singlet and triplet states of the mixtures in Eq. (11) contribute to the widths. If we use the perturbative 
wave functions, then the relative phase of the triplet contribution with respect to the singlet contribution is 7r/2. 
The widths in this case are 

/ \ 2 / \ 2 

Vfi= ^^L { V— -2 <2 yl .n T . < 2 50 + 7 (u) 

±7 




The resulting radiative widths are shown in Table III. 
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r 7 (keV) Model Expt 



D* a D s 


1.91 


< 1.9 x 10 3 


D s0 (2317) -> D s * 


4.92 




D s i(2460) -> D s 


12.8 


BR =0.18 ±0.04 


D s i(2460) ^D s * 


15.5 


BR< 0.08 


D s i(2460) -> D s0 (2317) 


5.74 




D^i(2536) D s 


54.5 




D' al (2536) ^ D* s 


8.90 


possibly seen 


D' sl (2536) -> D s0 (2317) 


2.36 




D a 2 (2575) ^ d; 


44.1 




D s *(2637) -> D s0 (2317) 


6.76 




D s *(2637) D s i(2460) 


2.8 




D s *(2637) D^(2536) 


0.24 




£> s *(2637) -> D s2 (2573) 


0.35 




£> s (2485) -> D s i(2460) 


0.01 





TABLE III: The the radiative decays of the D a mesons are shown. These widths are computed using the mass values 
obtained directly from our calculation. This includes the n = 2 pseudoscalar and vector states, the latter of which has 
recently been observed with a higher mass [13]. The widths are from [13]. 

5. CONCLUSIONS 

We have shown that a potential model consisting of the relativistic kinetic energy, a linear long-range confining 
potential together with its v 2 /c 2 relativistic corrections, and the full v 2 /c 2 plus one-loop QCD corrected short 
distance potential is capable of providing extremely good fits to the spectra of the D s states by treating them 
as states of the cs system. We find that in this perturbative treatment the long-range potential must be entirely 
due to scalar exchange. 

The single photon widths can be obtained from the variational wave functions, but, apart from some branching 
ratio measurements, there are relatively little data available. Our theoretical results are comparable to those 
given in Refs.[5] and [7] allowing for the fact that both of these references use a substantially higher strange 
quark mass (419 MeV and 480 MeV, respectively). In every case, efforts to model these states will be greatly 
improved by the availability of additional data. 
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APPENDIX A: DETAILS OF MIXING 

The mixing of the 3 Pi and 1 P\ states is obtained by diagonalizing the 2x2 matrix 

where E3 is the 3 P\ energy, E\ is the 1 P\ energy and V31 is the mixing matrix element. In perturbation theory, 
this is relatively simple since all of these matrix elements can be calculated using the unperturbed wave functions 
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that are all the same. The energy eigenvalues are 

E± = \ (E 3 + E 1 ) ± ly/(E 3 - E 1 r + W 3 \ , (A2) 

and we fit the D^(2536) and £> s i(2460) to E + and Note that as V 31 -> 0, E + -> £ 3 and E_ -> El To 
define the mixing angles in terms of known parameters, we assume that the eigenvector tf> + corresponding to E + 
behaves as 

^ V ^°(l)- (A3) 



With this assumption, is 



By writing ip+ as 



with 



(E + -E l Y + V 3 \ 



tan(fl) = - Vsi , (A6) 

£,+ — £j\ 

we arrive at the decomposition Eq. (11). It is possible to obtain an estimate of the mixing angle by using the 
branching ratios of the £> s i(2460) — > D s j and £> s i(2460) — > D* 7, whose ratio gives 



r(£ sl (2460) ^ D* s l)uJ 3 
T(D sl (2A60)^D sl ) u* 3 



tan 2 (6»), (A7) 



where w is the momentum of the photon in the D s transition and u>* is the corresponding photon momentum 
in the D* transition. Using the published branching ratio information, with favorable assumptions, the data are 
consistent with 9 = ±50°. Our calculation gives 9 = 59.7°. 
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